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Monotonicity of Dissipative Flow Networks Renders Robust Maximum Profit 
Problem Tractable: General Analysis and Application to Natural Gas Flows 


Marc Vuffray, Sidhant Misra, Michael Chertkov 


Abstract —We consider general, steady, balanced flows of a 
commodity over a network where an instance of the network 
flow is characterized by edge flows and nodal potentials. Edge 
flows in and out of a node are assumed to be conserved, thus 
representing standard network flow relations. The remaining 
freedom in the flow distribution over the network is constrained 
by potentials so that the difference of potentials at the head 
and the tail of an edge is expressed as a nonlinear function 
of the edge flow. We consider networks with nodes divided 
into three categories: sources that inject flows into the network 
for a certain cost, terminals which buy the flow at a fixed 
price and “internal” customers each withdrawing an uncertain 
amount of flow, which has a priority and thus it is not priced. 
Our aim is to operate the network such that the profit, l.e. 
amount of flow sold to terminals minus cost of injection, is 
maximized, while maintaining the potentials within prescribed 
bounds. We also require that the operating point is robust 
with respect to the uncertainty of customers’ withdrawals. In 
this setting we prove that potentials are monotonlc functions 
of the withdrawals. This observation enables us to replace 
in the maximum profit optimization infinitely many nodal 
constraints, each representing a particular value of withdrawal 
uucertaluty, by only two constraints representing the cases 
where all nodes with uncertainty consume their minimum and 
maximum amounts respectively. We illustrate this general result 
on example of the natural gas transmission network. In this 
enabling example gas withdrawals by consumers are assumed 
uncertain, the potentials are gas pressures squared, the potential 
drop functions are bilinear in the flow and its Intensity with an 
added tunable factor representing compression. 

1. Introduction 

The maximum profit, or alternatively minimum loss, net¬ 
work flow problem aims to maximize monetary benefit by 
delivering maximum amount of flow (of a commodity) from 
sources to terminals. The setting is general, and as such it 
applies to natural gas networks [1], [2], [3], our enabling 
example, but also to electric circuits [4], [5] and traffic flows 
[6], [7], [8], 

Formally, these problems are constructed by extending 
the standard network flow setting, see e.g. [9], [10] and 
references there in, with additional physical constraints, 
introducing nodal potentials and relating the potential drop 
along an edge of the network to a function of the flow. Thus, 
in the case of the gas flows, the potentials are pressures 
squared and the pressure square drop is a bilinear function 
of the flow and the flow amplitude with added term related 
to compression [1], [2], [3]. 
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In this manuscript we focus on discussing a robust version 
of the maximum profit network flow problem, more accu¬ 
rately “adjustable robust optimization” problem following the 
terminology commonly accepted in the literature on robust 
optimization [11], [12], [13], [14]. This means that in the 
robust optimization network flow model considered there are 
three different types of variables: the uncertain variables, the 
non-adjustable variables and the adjustable variables. The 
uncertain variables express information that is not certain, 
i.e. available for the optimization decision only in the form 
of allowed range. The non-adjustable variables represent the 
“here and now” decision in the system. Their values should 
be feasible for any realization of the uncertain variables from 
the allowed range. Finally the adjustable variables represents 
the “wait and see” decisions. Their values are adaptable to 
a particular values of the uncertain parameters. 

The adjustable robust optimization is composed of a test of 
robust feasibility and an optimization procedure. A value of 
the non-adjustable variables is said to be robust feasible if for 
any acceptable configuration of the uncertain variables there 
exists feasible values of the adjustable variables. Then the 
optimization procedure consists in finding a robust feasible 
configuration of the non-adjustable variables such that the 
objective function, maximum profit or minimum loss, is 
minimized. 

The difficulty in solving robust optimization problems 
arises primarily due to the robust feasibility constraint, which 
is in essence an intersection of infinitely (and possibly 
uncountably) many constraints, one corresponding to each 
allowed value of the uncertain parameters. This results 
in the so-called semi-infinite program [15]. In the robust 
optimization literature, the ways to handle these constraints 
can be classified into three different categories. First, when 
the constraints and the uncertainty set have special structure, 
e.g., linear constraints and ellipsoidal uncertainty set, it is 
possible to use duality theory to represent the infinitely many 
constraints with one single dual feasibility constraint [16], 
[17], [18]. Also included in this category is approximations 
and relaxations of more complicated uncertainly sets and/or 
constraints with simpler ones that are amenable to duality 
theory. The second category/approach is similar to the so- 
called “scenario based” approach, where a (possibly ran¬ 
dom) sampling of the uncertain parameters is performed, 
and the feasibility constraint corresponding to each sampled 
parameter is included in the optimization formulation [19], 
[20]. The quality of the solution thus obtained depends on 
the number of samples used and also on how the samples 
were chosen. The third case, which is the approach taken in 



this manuscript, is when one can analytically or numerically 
identify the “extreme-cases”, i.e, find the subset of values 
of the uncertain parameters that can violate the feasibility 
constraints. When this subset is finite, or has a finite repre¬ 
sentation, the robust feasibility constraint again reduces to a 
finite number of standard constraints. Examples where this 
strategy is used are scarce. (See [20, pp. 388] for a brief 
discussion on the topic.) 

Our main result, stated in Theorem[T] is that the adjustable 
robust maximum profit problem is tractable, in the sense 
that instead of keeping infinitely many conditions, associated 
with all possible realizations of the uncertain variables, it is 
sufficient to only account for two extremal conditions corre¬ 
spondent to every uncertain variable (customer withdrawal) 
greedily maximize/minimize their values. We implement a 
multi-stage strategy to prove the results. First, we prove in 
Theorem[^existence and uniqueness of the optimal solutions 
for the adjustable variables given that the other parameters in 
the problem are fixed. Then, we prove in Theorem]^ that the 
adjustable variables are monotonic functions of the uncertain 
variables. Finally, we combine all these result to prove the 
main tractability statement. Schematic outline of our proof 
strategy is shown in Fig. [T] 
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Fig. 1. Schematic diagram relating different statements of the manuscript. 


The manuscript is organized as follows. In Section]^ we 
briefly introduce the concept of dissipative network flows 
and define the problem of the adjustable robust maximum 
profit/minimum loss problem. In this Section we also state 
our main theorem/statement. In Section III we prove exis¬ 


tence and uniqueness of a configuration of the adjustable 
variables given that the other parameters are fixed. The 
proof is based on the energy function method which consist 


in rewriting a set of equations as the result of a convex 
minimization problem. In Section IV we prove that the ad¬ 
justable variables are monotonic with the respect to uncertain 
variables. This result is general - it applies to the case 
of unconstrained network flows. Finally, in Section [V] we 
discuss application of our results to robust optimization in 
natural gas transportation networks. 


II. Formulation and the Main Result 
A. Dissipative Flow Networks 

Fet G = {V, E) be a connected directed graph, where 
V and E C V X V denotes the set of nodes and directed 
edges respectively. Our notational convention is that for 
any pair of , i and j, connected by an edge, both {i,j) 
and (j, i) belong to E. Function g ; V —)■ K associates 
production/consumption qi to the node i G V. We assume 
that the total production/consumption is balanced 

Y.q^ = 0■ ( 1 ) 

iev 

Flow (/) : F —>■ R is a function mapping an edge {i,j) G 
E to flow intensity (j)ij (which we also call flow, or edge 
flow, when not confusing). The flow is skew-symmetric with 
respect to the edge inversion, i.e. for all (z,j) S E 


^ij — ^ji- ( 2 ) 

The flow is conserved locally, i.e. at any i G V equation 

(t>ji + qi = o, ( 3 ) 

j^di 

where we denote the set of neighbors of i G V hy di := 
{j G V \ {j,i) G E}. Note that in the context of electric 
circuits, the flow conservation equation (|^ is also called the 
first Kirchhoff’s Faw. 

Flow over an edge (i,j) G £ is related to the potentials at 
end-nodes of the edges through the following potential drop 
equation 

Tj T* = fij , (4) 

where the functions, coined dissipation functions, are 

assumed to be continuous and monotonically decreasing, 
thus leading to decrease/dissipation of potential along the 
direction of the flow. The dissipation functions satisfy the 
following symmetry relation with respect to an edge inver¬ 
sion 

fij (x) = -fp (-x), (5) 

and they also possess an operational freedom, i.e. the dis¬ 
sipation functions can take any form from their admissible 
set, f G E depending on the problem considered. 

Note that existence of a potential function that satis¬ 
fies Eq. 0 allows a graph cycle interpretation. Fet 7 = 
{ii,..., in} G V" be a directed path in G such that ik ^ ii 
for fc, I G {1, ...n — 1}. If ii = in we say that 7 is a directed 
cycle, otherwise we called 7 a non-intersecting directed path. 
For any directed cycle 7 in G we sum the potential drop 






















Eq. 0 along the directed edges of 7 to get the following 
expression 

n— 1 

fikik+l {4>ikik+l) ~ ( 6 ) 

fc=l 

Eqs. 0 can be interpreted as a nonlinear generalization of 
what is called second Kirchoff’s law in the circuit theory. 
In fact, conditions expressing the fact that the flows satisfy 
Eqs. 0 for any cycle basis of the graph is fully equivalent 
to the set of conditions 0 over all edges of the graph. 

We call the network G = {V,E), with flows, potentials 
and the dissipative functions, /, defined on it, the dissipative 
flow network. 


B. Robust Maximum Profit/Minimum Loss Optimization 


We consider a dissipative flow network where the set of 
nodes is partitioned into three non-empty subsets 1/ = ^ U 
T U R. The set S contains nodes i G S that are sources 
injecting flow into the network which costs gi (qi), where gi 
is a pre-set function. Nodes i G T are terminals where the 
flows are withdrawn leading to the following payment per 
node, hi (q^), where hi is non-decreasing function. Nodes i G 
R are internal customers each consuming a fixed amount of 


flow that is uncertain, i.e. it lies in the set q^,qi . We assume 
that the internal consumption is mandatory, and as such it is 
not priced. We aim to maximize the revenue, defined as a 
difference between the profit brought by selling the flow at 
the terminals minus the cost of injection. Alternatively we 
can state the problems as minimization of expenses, defined 
as the cost of injection minus the profit at the terminal: 


c{qs,qT) -^hi{qi) , (7) 

iGS iGT 


where qs and q^ are shorthand notations for production at 
the sources and the terminals respectively. 

The maximum profit (minimum loss) problem is an un¬ 
certain program that minimize the cost 0 with respect to 
injections, potentials, flows and dissipation functions from 


their admissible set R 


minimize 

c{qs,qT) 


subject to 

realization 

Hi — Hi 

Vi G S, 


^ ^ 4^ji Qi — 0 
j^di 

ViGV 


7Tj — 7Ti = —fij 

^* = 0 
iev 

V (i, j) e E 


TTi < TTi < Tfi 

Vi G C, 


where the regular customer withdrawals could be any, but 
fixed, value ^realization g q., g ■ 


The optimization 0 is deterministic as the uncertain 
parameters are fixed. Below we will turn to its robust equiv¬ 
alent, thus aiming at taking optimization decision leading to 
a solution which is feasible for any realization of uncertainty 


while also achieving losses which are first maximized (ad- 
versarially) over the internal (uncertain) degrees of freedom 
and then minimized over the external (operational) degrees 
of freedom. Therefore, to arrive at a plausible robust opti¬ 
mization formulation one, first, need to split all the variables 
into the sets of the so-called operational (non-adjustable) 
and adjustable variables [12]. Operational variables represent 
decisions taken regardless of a particular configuration of the 
uncertainty. Adjustable variables, according to the name, are 
flexible, i.e. they adjust themselves to the uncertainty. 

We consider the case where the operational variables are 
the production values at the source qs := {qi}i^S’ 
potential values at the terminals := {TTi}^^j, and the 
dissipation functions f G R. The adjustable variables are 
the production values at the terminals qr := {qi}i^T’ 
potential values at the source tts ■= {TTij^gg, the potential 
values at the consumers Tr^ := {TTi}^^j^ and the flows, 
(j). The production values at the consumers are denoted by 
qR '■= ^nd the cartesian product of their uncertainty 

set is denoted as Q = g. x ••• x g. ,g,- , . We 

refer the reader to the illustration of the natural gas network 
in Section to justify this particular choice of adjustable 
and non-adjustable variables. 

With this choice of the operational variables, an adjustable 
robust counterpart of the deterministic minimum cost prob- 


lem 0 becomes 


minimize 

qs,'^T,f,x 

X 


such that 

VqR G Q, 3qT,Trs,TTR,(l> 


subject to 

c{qs,qT) < X 



fji + gi = 0 

jedi 

Vi G C 


'Xj — TTi = — fij ifij) 

y]g* = o 

i&V 

V(i, j) G E 


Tfi < TTi < 7fi 

Vi G V. 


(9) 


Considered directly, Eq. 0 is intractable as containing 
an infinite number of constraints, each associated with a 
particular qs G Q. However, and in spite of the grim naive 
assessment, Eq. 0 allows tractable re-formulation stated in 
the following main theorem/result of the manuscript: 


Theorem 1 (Tractability of the Robust Minimum 
Loss/Maximum Profit Problem). There exists a function 
q^ : kI-R|+I'S'I+|7’| x R —>■ and a function 

~RUS . ]^|i?|+|s|-r|T| X R ^ r\Rus\ 

adjustable robust counterpart 0 reads 
minimize c{qs,qT {qR,qs,'^T, f)) 

subject to TTj < •^) Vi G i? U S' 

(9K>9s,7rT,/) < 7fi ViGSUS 
ZTi S < ^i Vi G T, 


where g„ := 

—rt 



ieR 


and q^ := {gj^g^ 


are respectively 










the minimum and maximum withdrawal values at the con¬ 
sumer nodes. 


The functions Tf" and mentioned in Theorem are 

derived explicitly, moreover stated in terms of the strictly 
convex optimizations, in Section |ng In words Theorem [T] 
tells us that only the two extreme realizations of have to 
be considered in order to solve the semi-infinite optimization 


©■ 

Our strategy to prove the Theorem is the following: We 
first prove that there exists a unique solution of the adjustable 
variables qx, tts, t^r and (j) given a configuration g/j, qs, 
ttt and /. This states the existence of the functions if" 
and . Then, we show that the solution qx, tts, ttr are 
monotonic functions of qu, that enables us to reduce number 
of critical constraints to only two corresponding to the 
extreme configurations of the uncertainty. This establishes 
that f" and are monotonic functions of qR. 

The proof of existence and uniqueness of the solution is 
done in Section III using the energy function method. The 
monotonicity properties are proved in Section IV and only 
require some properties of the flow network. 


III. Uniqueness oe Dissipative Flow Solutions: 
THE Energy Function Method 

In this Section we show the existence and uniqueness of 
a configuration of adjustable variables, given a configuration 
of the non-adjustable variables. The key idea here is to relate 
solutions of the dissipative flow network equations ([T]i, Q 
and to the extremum of some convex “energy” function. 

Theorem 2 (Uniqueness and Existence of the Dissipative 
Flow Solutions). Let G = (U, E) be a dissipative flow 
network. Given qn, qs, ttx ond f, there exists a unique 
solution qx, tts, ttr and (j) which satisfies the balanced 
production equation Q, the flow conservation equation Q 
and the potential drop equation 0- 

Proof. Let us first invert the potential drop Eq. 0 to express 
the flow with respect to the potential 

fij ~ '^j) = 4>ij- 

This operation is possible because is an increasing 
function. Then, we rewrite the flow conservation equations 
only in terms of the potential 

+ 9* = 0 . 

j&di 

Therefore if one finds ttr and ttr that satisfy the flow 
conservation equation for i G S U R, one reconstructs the 
production qx using flow conservation for i G T. It is easy 
to see that the balanced production Eq. 0 is automatically 
satisfied 

iGV i^V jGdi 

= - 

ihpeE 

= 0 , 


where in the last line we use the symmetry relation 0 of 
the drop function. 

We use the energy function method to prove that the po¬ 
tentials TTS and ttr are uniquely determined. Let us introduce 
the set of oriented edges 

O ■■= G E\i <j}. 


The graph T = (U, O) is an orientation of G i.e. {i,j) G T 
if and only if {i,j) G G and {j,i) ^ T. 

Consider the following energy function 


£ {tTS,TTr\ qR,qS,TTT,f) = ^ {tTi - TTj)-'^TT,qi, 

(jj)GO 

( 10 ) 

where Eij is a primitive of f~f. Note that the primitive exists 
because is continuous as fij is continuous. It is easy to 
see that a minimum of £ with respect to ttr and ttr satisfies 
the flow conservation equation at z G StJR. Indeed by taking 
the derivative, one arrives at 


d£_ 

Ott, 


j^di 


We now have to prove that the extremum is unique and that 
the extremum is actually the minimum. 

Let us define for every edge (z,j) G O, variables 


•— TTi TTj . 


It is straightforward to see that the function 

i^ij) ) 

(ij)eo 

is a strictly convex function of A^. Note that is 

increasing because fj is increasing. This implies that T'ij 
is also a strictly convex function of A^. Now we have to 
prove that it the function is also strictly convex in ttr and 
ttr. 

First we express the relation between A^ and tt^ in 
a matrix form using the incidence matrix M of E. The 
incidence matrix M is a \V\ x \0\ matrix with entries 
= 5kj — Ski where 6 is the Kronecker’s symbol. 
The relation between A = {Aij}^^^Q and tt reads 

A = —M'^tt. 


We can also explicitly separate contributions of the potential 
TTR and TTfj from ttx in this equation by introducing the 
reduced incidence matrices My\x and Mx. The matrices 
My\x and Mx are equal to M with the lines, corresponding 
respectively to to z S T and i G V\T respectively, removed. 
One Anally obtains the following equation 

Observe that the matrix Myyj, is full rank. To see this, note 
that the product My\xMy\x corresponds to the Laplacian of 
the graph L = MAft with the lines z and j, corresponding to 
i,j G T, removed. Due to the Kirchoff’s matrix-tree theorem 




[4] we know that the determinant of any co-factor of L is 
non-zero on a connected graph, which implies that 

rank(MT^^) =|y\r|. 

This proves that S is strictly convex in tt^ and □ 

Theorem [^enables us to eliminate the adjustable variables 
from the adjustable robust counterpart We can now 
rewrite our optimization problem as follows 


minimize x 

9S.7rT,/,a: 

such that V( 7 ij S Q, 

subject to c {qr, ^ {qr, qs, ttt, /)) < x 

7L^ {qR,qs,'^T,f) <^i yiGR,S 

TTj < TTi < Tfi Vt G T. 

( 11 ) 

Function gs, /) which outputs the potential 

TTi for i G RU S was explicitly constructed in the result of 
the strictly convex optimization 

(te, qs, TTr, /) := arg min £ (tts, ttu \ q^, qs, t^t, f), 

( 12 ) 

where £ is defined by ([Tgi. Once the potentials are found 
using Eq. one can easily reconstruct the productions 

(qR, qs, ttt, f) = Yl ~ 

j^di 


IV. Monotonicity Properties of the Dissipative 
Flow Networks 


The existence and uniqueness properties from Section III 
allowed us to simplify the robust minimum loss optimiza¬ 
tion problem. However, it still remains in the form of an 
intractable semi-infinite program given by Eq. O- In this 
Section we show that only two scenarios for qR has to be 
considered in order to solve Eq. GB- To achieve this goal 
we show that the potentials given by Eq. and the 

productions given by Eq. G3 are monotonic functions 

of qR- 


A. Flow Networks 

The following theorem applies to general network flows, 
i.e. flow networks which are not necessary dissipative. We 
only require here that the flows satisfy the flow conservation 
Eq. This very general result will then become a starting 
point to prove monotonicity properties of the dissipative flow 
networks. 


Theorem 3 (Aquarius Theorem). Let G = (E, E) be a 
flow network and let (j) and (jf be flows that satisfy the 
flow conservation Eq. for the productions q and q* 
respectively. Let T <Z V be a subset of V. If qi > q* for 
all i G V \ T, then for every node u G V \ T there exists a 
non-intersecting path {A, ■ • •, in} such that ii G T, in = u 
and Moreover if qu > qu the inequality is 

strict i.e. 

Proof. We construct the path by induction. Choose u G V\T 
and assume that qu > ql^. The proof in the case where qu = 


qu is identical. We define a sequence of subsets of nodes 
Bj, CV and Bi in the following way 

Bi := {m} 

Bk+i ■= {* S E \ Ak I GdiFBk s.t. 4)*^ > fij] . 

An example of the sets Ak and Bk is shown in Fig. 



Fig. 2. Constmction of the ensembles Af^. Note that B]^ = A^ \ Ak_i. 
The edges for which (p*- > (pij with i E and j E B]^ are shown 

with a solid line and a dark arrow. 

Let n be the first k for which BkDT First, we prove 
by induction that the sets Bk for k < n aie non-empty. The 
set Bi is non-empty by construction. Suppose that the sets 
Bk are non-empty. Define the set of edges connecting two 
vertices in Ak by 

Ek := {{i,3) e E I z,j S Ak} , 

and the set of edges with a starting point in E \ and an 
endpoint in Ak 

^Ak := {{j,i) G E \ j ^ Akandi G Ak}. 

By summing the continuity Eq. 0 for the flow f over the 
vertices in Ak one obtains 

0 = H ) 

ieAfc y jGdi J 

= 5Z 51 5Z 

ieAk ieAi, jedi 

= qi+ 51 + 4>jt) + 51 

= 5Z ^ 5Z ^p’ 

where the skew-symmetry of the flow 0 between the two 
last lines has been used. Summing up the same continuity 
relations for the flow </>* leads to the following inequality 

51 = - 51 * < - 51= 51 ^*p- 

U,i)etA,, teAk ieAk {j,i)etA,, 

The above inequality implies that there exists {v, w) e ~^Ak 
such that fuw > 4>vw It remains to be checked that the node 






w G Ak is an element of Bk = Ak \ A^-i- By construction 
if w S Ak-i then v should be contained in Ak which 
is a contradiction. Thus Bk+i is non-empty. We have just 
proved by induction that B^. for fc < n are non-empty and 
moreover that Ak is an increasing sequence of sets. The last 
observation guarantees that n is a finite number. 

We now focus on proving that there exists a path 
from some t G T to u such that > 

By induction it is now easy to prove that for every 
node i G Bk there exists a path of length k from i to u 
such that > 'Piiii+i- Since it was already shown that 

there exists n S N such that Bk D T 7 ^ 0 the proof is 
concluded. □ 

Theorem has an intuitive explanation, clarifying the 
name chosen for the theorem. Picture a flow of water through 
a network of aqueducts. If water leaks appear at some points 
within the network the system becomes unbalanced. There¬ 
fore an aqueduct inspector (“Aquarius” in Latin) needs to 
compensate for the leaks pouring more water at other nodes. 
This extra water added at the injection node will flow along 
some path that has to end up at a leak. Otherwise the added 
water accumulates and the network remains unbalanced. 

B. Dissipative Flow Networks 

Theorem [^applied to the dissipative flow networks results 
in the following statement. 

Corollary 4 (Monotonicity of Potentials with Consump¬ 
tions). Let and (</>*, tt*) be solutions of the continuity 

Eq. 0 and the potential loss Eq. 0 for the productions q 
and q* respectively and the same dissipation function f. If 
ttt > Ttt for all t G T C V and qi > q* for all i G V \ T 
then TTu > TT*/or every node u G V\T. Moreover if q^ > g* 
the inequality is strict i.e. 7r„ > tt*. 

Proof Consider u G V\T. One evaluates the potential drop 
Eq. Q as explained by Theorem to arrive at the inequality 

71 

k^l 

71 

k=l 

where we have just used (in between the two last lines) that 
fij is an increasing function. □ 

Corollary allows us to make the following statement 
about relations between productions and consumptions. 

Corollary 5 (Monotonicity of Productions with Consump¬ 
tions). Let (/, tt) and ((/)*,7r*) be solutions of the continuity 
Eq. 0 and the potential loss Eq. 0 for the productions q 
and q* respectively and the same dissipation function f. If 
TTt = TT^ for all t G T G V and qi > q* for all i G V \ T 
then qt < q^ for all t G T. 


Proof. Choose u GT and assume that qu > Qu- Corollary]^ 
applied to u implies that tt^ > tt* which is a contradiction. 

□ 

Corollary directly implies that the maximum (and 
respectively minimum) of the potential iTi {qji, qs,TTT, f) 
given by Eq. OH is achieved at a maximum (respectively 
minimum) of the customers production qn 

max TTj (qij,gs,7rT,/) = tt, (g^, qs, ttt,/) 
ir&Q 

min TT, (gij,gs, 7 rr,/) = tt* gs, ttt,/) . (14) 

With regards to the Corollary]^ it tells us that the production 
Qi {qR, qs,TtT, f) given by Eq. OH achieves its maximum 
(respectively minimum) at g^ (respectively at g^). This 
implies that the maximum of our objective function is 
achieved when all the regular customers are consuming at 
their maximum value 

{qR, qs, TTt, /)) = X {qi {qR, qs, t^t, /)) • 

(15) 

Here in Eq. ( fTS] ) we use the property that ft, is a non¬ 
decreasing function of g^. 

Eqs. ( [T4| and ( [T5| l show that only the two extreme 
scenarios of the regular customers consumption has to be 
considered to guarantee feasibility of all other consumption 
configurations from the uncertain intervals. 

V. Illustration with Natural Gas Networks 

A natural gas network is a system of interconnected 
pipelines delivering natural gas from producers to consumers. 
In the normal operational regime gas flows in the transmis¬ 
sion (high pressure level) pipes are turbulent. Consider a 
pipeline of length L. As a pipeline length is much longer 
that its cross-section, the system can be modeled as one¬ 
dimensional, parameterized by position along the flow, x G 
[0, L], with cross-section effects averaged out. Then, the state 
of the gas flow at a time t is characterized by its pressure 
p {x, t) and mass flow / (x, t) along the pipe (both averaged 
over the cross-section). The two characteristics are related 
to each other via the following set of partial differential 
equations [1], [3], [21], [22] 

dt(j){x,f) + da,p{xA) = 

2p [x, t) 

dtp{x,t) + cldx(j){x,t) = 0 , (16) 

where Cs is the speed of sound in the gas and a is a constant 
that depends on the type of gas, size of the cross-section, 
roughness of the pipe surface and also on the Re-number, 
characterizing the level of turbulence. The term K{x,t) in 
Eqs. ( [Thl l accounts for a compression added at a compression 
station to compensate for pressure drop. The first equation 
in Eqs. ( [Thl l is a phenomenological equation that quantify 
the loss of momentum due to turbulent friction. The second 
equation in Eqs. ([T^ enforces conservation of the fluid mass 
along the pipe. Compression can be modeled as 

K {x, t) = bS {x — Xc) , 


(17) 



where 6 (•) is the Dirac’s delta function; Xc is the compressor 
station position along the pipeline and b is an additive 
compression factor of the station. 

The equations are dynamic, however for the purpose of 
planning the gas flow budget on the scale of a day, the 
dynamics in Eqs. can be ignored, thus setting dtip = 
dtp = 0. In this case, after straightforward spatial integration 
along the pipe, Eq. transforms into 

p{L,tf -p{0,tf = -^(j)\(j)\+b 

(j){x,t) = (j), (18) 

relating pressures at the ends of the pipe to the amount of 
flow and the value of compression acquired along the pipe. 
Eq. ( [T8 ] i tells us that the flow in one pipe is constant and 
that it is driven by a difference of pressure squared at the 
endpoints of the pipeline. This static representation enables 
us to model the gas network in the steady-state regime as 
a dissipative network with a potential equal to the pressure 
squared tt = p^. The dissipation function over a pipe (z, j) 
is nonlinear; 

L'' Oi '' 

fij ii’ij) ~ ^ I'/’ul + (19) 

The three types of nodes that we have introduced above 
in the maximum profit problem (sources, internal customers 
and terminal) map into nodes of the gas network as fol¬ 
lows. The sources that inject gas into the network are 
gas producers, e.g. gas processing plants. Liquid Natural 
Gas (LNG) terminals and storage injecting gas into the 
system. The internal customers with uncertain demand are 
consumers with existing contracts for gas delivery, such as 
Local Distribution Companies (LDC) and electric gas-bred 
plants. The terminals are opportunistic customers ready to 
buy whatever amount of gas which can be made available 
(e.g. LNG terminals and storage reservoirs working in the 
regime of gas accumulation). The uncertainty on the side of 
the internal customers accounts for exogenous changes such 
as those related to LDC consumers’ heating requirements, 
and fluctuations of gas consumption at the gas-bred plans 
due to uncertainty on the electric grid side (e.g. the renewable 
generation). 

The choice of operation variables is not unique and de¬ 
pends on the regime of operation, type of the system, country, 
etc. Lor example in US, typical large-scale producers would 
maintain constant injection/flow, thus allowing changes in 
the pressure, while LDC (and related city-gates) on the 
contrary would withdraw constant flow allowing the pressure 
to meander. 

VI. Path Lorward 

In this manuscript we proved that the robust maximum 
profit problem over nonlinear dissipative network flow prob¬ 
lem is tractable. The strategy that we have employed in the 
proof is based on the search for an explicit formulation of 
the special limiting scenarios such that feasibility of solution 
for the special scenarios guarantees feasibility for all other 
scenarios from the uncertainty range. We proved that in 


the general case of the static dissipative network flow it is 
sufficient to maintain feasibility only for two scenarios. The 
essence of this major step in our proof strategy is related 
to the very strong monotonicity property in the space of 
solutions. 

We envision extending this work in the future along the 
following three directions: 

• We plan on moving from static formulation to dynamic 
and thus to analyze dynamic versions of respective 
robust optimizations from the perspective of scenario 
reduction discussed in this manuscript. There are more 
than one possible generalization strategies (for transition 
from static to dynamic). In particular, one may hope 
to And a dynamic (Lagrangian) generalization of the 
Aquarius principle, i.e. generalization of the Theorem 

In the context of natural gas application, we plan 
to introduce and analyze consumption robust version 
of the dynamic optimization (off-line control) problem 
discussed in [23] based on the dynamic Eqs. ( [T8) l. We 
will also attempt to develop dynamic version of our 
scenario-reduction technique suitable for control and 
dynamic optimization traffic network problems of the 
type discussed in [6], [7], [8]. 

• We would like to design an efficient numerical scheme 
to solve the tractable version of the robust maximum 
profit/minimum cost problem. As formulated in Theo¬ 
rem the problem is a bi-level optimization task. To 
advance this task we plan to utilize the energy function 
representation in order to formulate the entire problem 
as one minimization procedure in the spirit of [24]. 

• Linally, we plan to extend the results reported in the 
manuscript to more general types of uncertainty sets 
such as ellipsoids. The ellipsoid type of uncertainty 
set is more challenging than the one of the box kind 
considered in the manuscript. Even if the monotonicity 
properties guarantee existence of the extremal scenarios, 
deriving explicit form of the extremal scenarios in 
setting other than of the box type remains a challenge. 
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